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Abstract—The first results obtained from the aerosol-cloud re- 
trieval algorithm (developed at the University of Bremen) are pre- 
sented. The algorithm enables the observation of the regional char- 
acteristics of aerosol and cloud optical thickness both over land 
and ocean surfaces. The aerosol and cloud optical thickness over 
Western Europe is derived from the high-resolution SeaWiFS data 
for October 11, 2001 (11:30 UTC). The most probable value of the 
aerosol optical thickness was found to be equal approximately 0.25. 
The frequency distributions of the aerosol and cloud optical thick- 
ness are skewed and have long tails for larger optical thickness. It 
was found that retrieved values of the aerosol optical thickness at 
wavelengths 0.412 and 0.440 um are close to those measured by five 
ground-based instruments placed at different locations. The prob- 
lems related to the retrieval of the atmospheric optical thickness 
from space are discussed. 


Index Terms—Aerosols, atmospheric optics, clouds, radiative 
transfer, remote sensing. 


I. INTRODUCTION 


LOUDS AND aerosols are important constituents of the 
C atmosphere and in particular of the troposphere. In addi- 
tion to their unique role in the hydrological cycle of the earth—at- 
mosphere system, they act as both sources and sinks for atmo- 
spheric gases in both the stratosphere and troposphere as well as 
being of significance in determining the radiative balance within 
the atmosphere. In particular, cloud processes supply and de- 
plete water vapor. Our knowledge about the global behavior of 
clouds and aerosols is poor and not sufficient to assess the im- 
pact of anthropogenic activity on atmospheric chemistry and cli- 
mate change (e.g., see the Intergovernmental Panel on Climate 
Change (IPCC) [24]). 

The remote sensing of aerosol and cloud fields from the in- 
strumentation aboard space-based platforms provides a unique 
opportunity to improve our understanding of the global behavior 
of clouds and aerosols [6]. In particular, the optical thickness 7 
of terrestrial atmosphere as a function of the wavelength results 
from these measurements. In the visible spectral region, where 
the energy of the incident solar flux has a maximum, this pa- 
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rameter is dominant in determining the radiative energy budget 
of our planet [57]. 

The atmospheric optical thickness, excluding those spectral 
regions where atmospheric gases exhibit strong absorption, de- 
pends mostly on the column concentration of suspended parti- 
cles within the atmosphere: the latter comprising, e.g., dust and 
water-soluble aerosols, cloud droplets and crystals. As 7 also 
depends on the size, shape and chemical composition of parti- 
cles, the slope of the aerosol optical thickness as a function of 
the wavelength is an important characteristic, which yields the 
aerosol type. 

Several studies have been devoted to the derivation of the 
aerosol and cloud optical thickness from measurements of 
the top-of-atmosphere reflectance, performed by different 
radiometers and spectrometers on board a variety of satellite 
platforms [10], [13], [14], [22], [25]-[27], [31], [32], [41], 
[42], [49]. Aerosol and cloud optical thickness retrievals have 
mostly been performed by the different groups using different 
retrieval codes, designed for the retrieval of either cloud or 
aerosol parameters. However, clouds and aerosols exist in the 
terrestrial atmosphere simultaneously and also are mutually 
interdependent [6], [18], [44]. For instance, Breon et al. [6] 
speculate that there is probably an anthropogenic impact on the 
cloud microphysical and, therefore, radiative properties. This 
is, in particular, due to the influence of anthropogenic aerosols 
on cloud droplets spectra and concentrations. The increased 
droplet number concentration leads to smaller droplet sizes 
that make precipitation formation more difficult producing 
larger water contents [1], [17]. Clouds are not passive species 
in this interaction process. They also influence the aerosol 
microphysical properties [45]. 

In summary, it is, therefore, of importance to have combined 
measurements and retrieval algorithms, which are capable of the 
simultaneous derivation of both cloud and aerosol optical thick- 
ness (and other aerosol and cloud properties) both over ocean 
and land. Note that current aerosol operational algorithms are 
mostly restricted to the aerosol optical thickness retrieval over 
ocean. 

Algorithms, which are capable to derive the cloud optical 
thickness and separately the aerosol optical thickness for the 
same partly cloudy pixel currently do not exist. 

Our aerosol—cloud algorithm is also not capable to deal 
with partly cloudy pixels. However, we combined independent 
aerosol and cloud retrieval algorithms in one processing chain 
to have a possibility to derive simultaneously aerosol and cloud 
information in the scene under study (not in a single pixel) using 


0196-2892/04$20.00 © 2004 IEEE 


KOKHANOVSRY et al.: ATMOSPHERIC OPTICAL THICKNESS OVER WESTERN EUROPE 825 


a single retrieval procedure. The above-mentioned algorithms 
have been recently developed by the aerosol and cloud research 
group at the Bremen University for use with data from the 
Scanning Imaging Absorption Spectrometer for Atmospheric 
Chartography (SCIAMACHY) aboard the Environmental 
Satellite (ENVISAT) [7]. ENVISAT was successfully launched 
on March 1, 2002, and first data have already been delivered to 
the research community. 

The objective of this paper has been to demonstrate some 
of the capabilities and limitations of the algorithm developed 
for the studies of the atmospheric optical thickness. For this 
purpose, the well-proven data from the Sea-viewing Wide 
Field-of view Sensor (SeaWiFS) on Orbview-2, distributed by 
NASA Goddard Distributed Active Archive Center, have been 
used. SeaWiFS measures earthshine radiance in eight channels 
(0.402-0.422 um, 0.433-0.453 pm, 0.480-0.5 um, 0.500-0.520 
um, 0.545-0.565 jum, 0.660-0.680 um, 0.745-0.785 um, 
0.845-0.885 um). The swath width is 2801 km (58.3°) for the 
local area cover (LAC) product and 1502 km (45°) for the global 
area cover product (GAC). The spatial resolution is 1.1 km for 
the LAC and 4.5 km for the GAC. Note that the SeaWiFS was 
successfully applied to ocean color studies [15], [16], [58]. 
Stamnes et al. [53] introduced an algorithm for the simulta- 
neous retrieval of aerosol optical properties and chlorophyll 
concentrations over oceans using SeaWiFS imagery. 

A section of a global area cover data set of Western Europe 
taken on a cloudy day of October 11, 2001 was selected. 
It consists of 199 x 500 pixels with the size approximately 
4.5x 4.5 km? at the nadir observation. Thus, the results 
obtained are statistically significant. The correspondent file 
can be downloaded from the NASA Goddard Distributed 
Active Archive Center (file S2001284113048Lla_gac at 
http://daac.gsfc.nasa.gov/data/dataset/SEAWIFS/). 


II. AEROSOL OPTICAL THICKNESS AND ANGSTROM EXPONENT 


A detailed description of the aerosol subsection of the re- 
trieval algorithm is given elsewhere [56], and therefore, only a 
brief outline of the retrieval procedure is given below. 

The aerosol algorithm determines the aerosol reflectance 
R,(A) as the function of the wavelength by subtracting the 
contributions for the Rayleigh path radiance R,(A) and the 
contribution from the surface (land or ocean) reflectance R, (A) 
from the top-of atmosphere (TOA) reflectance Rypoa()). 
Namely, we have: Ra(A) = Rroa(A) — RA) — Rs (A), 
where Rs(A) = (Asurttit2)/(1 = Asur). Here, Ásurf is the 
spherical albedo of the surface, r is the spherical albedo of the 
atmospheric layer, and t1, t2 are transmission coefficients for 
light propagating from the sun to the surface (t1) and from the 
surface to satellite (t2). Functions tı and ty include both direct 
and diffused light transmission [33]. 

Lookup tables for the aerosol reflectance yield the spectral 
aerosol optical thickness Tą over oceans (A < 0.86 um) and 
vegetation (A < 0.67 um). Constrains are used to minimize 
the deviation of the aerosol optical thickness from a monotonic 
spectral behavior. We also find the parameters of the Angstr6m 
exponential law 7,(A) = GA~%, minimizing the functional y = 


(Ta(A) — BA~*)?. 


The algorithm takes into account the digital elevation model 
to determine the actual surface pressure for the determination of 
the Rayleigh path reflectance of the observed target [56]. This 
will be combined in the future with the retrievals using O2 and 
O4 absorption bands, which are, for example, available from 
SCIAMACHY [7] data. 

The crucial point in any aerosol retrieval scheme is the account 
for the surface reflection [29], [30], [58]. This is an important 
issue for aerosol retrievals both over land and ocean. In particular, 
Kaufman et al. [29], [30] proposed to use the correlation between 
spectral reflectances in visible and near-infrared (e.g., at 2.1 um) 
to solve the problem of the aerosol optical thickness determina- 
tion over land. This approach is routinely used by the MODIS 
team [8], [32]. This technique, however, cannot be applied to 
the SeaWiFs instrument due to the lack of the 2.1-jsm channel. 

Although it is of a general interest to be able to decode the 
SeaWiFS TOA reflectance in respect to the aerosol optical thick- 
ness both over land and ocean. 

Our correction for the land surface reflectance is based on the 
linear mixing of an apparent spectral surface reflectance from 
a “vegetation” spectrum and a “bare soil” spectrum tuned by 
the normalized differential vegetation index (NDVI). The dif- 
ferent surface spectra used were obtained during the LACE-98 
experiment [56]. 

A measured clear sky backscattered spectrum of the ocean 
is used to describe the light reflectance from water surfaces 
[56]. In particular, similarly to the case over land, we present 
the ocean surface reflectance as a linear mixture of experimen- 
tally measured clear water and coastal water spectra tuned by 
the normalized differential pigment index (NDPI), defined as 
the the differences of TOA reflection functions at wavelengths 
0.443 and 0.555 um divided by the value of the reflection func- 
tion at 0.49 um. The influence of whitecaps is neglected. Due to 
a comparatively small aerosol loading in the scene under study, 
the values of r and also the diffused parts of functions tı and 
t2 are small. Therefore, we keep only direct components of the 
solar flux in transmissions tı and tə. This simplifies the re- 
trieval. Note that the glint is avoided by the tilt of the instrument. 

One can find the flowchart of the procedure in Fig. 1. In par- 
ticular, note that the cloud free and totally cloud covered pixels 
are selected and then the optical thickness related to this par- 
ticular pixel is retrieved. The threshold value of the top-of-at- 
mosphere reflectance function R for the aerosol-cloud pixels 
discrimination is set to be equal to 0.2 at the spectral range 
0.845-0.885 jum (the eighth channel of the SeaWiFS). Values of 
R less than 0.2 correspond to cloud-free pixels in the algorithm. 
This is a rather crude estimation as compared to comprehensive 
state-of-art cloud screening algorithms (e.g., we do not account 
for the influence of the observation geometry on the threshold 
value of the top-of-atmosphere reflectance function, etc.). This 
could lead to the classification of thin cirrus clouds as aerosols 
and, therefore, to produce wrong retrieval results. However, it 
should be remembered that the SeaWiFS provides us with only 
limited spectral information in comparison with more advanced 
sensors (e.g., MODIS). Therefore, it should be underlined that 
results for the aerosol optical thickness (especially, for larger 
values of Ta) could be biased due to the unidentified presence 
of clouds in the field of view of the instrument. 
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Fig. 1. Flowchart of the retrieval procedure for the aerosol optical thickness 
determination. 


A large number of aerosol particles in the terrestrial atmos- 
phere (mainly of mineral origin and sea salt) are irregularly 
shaped. The use of the Mie theory, which is the basis of the most 
modern operational aerosol retrieval algorithms, therefore, may 
result in an inaccurate estimation of the aerosol amount. One im- 
portant feature of the algorithm developed is the precalculation 
of lookup tables using the experimentally determined aerosol 
phase functions for typical atmospheric conditions [56]. These 
experimentally determined phase functions include the features 
of light scattering by both spherical and nonspherical (e.g., pro- 
nounced side light scattering) aerosol particles. The semiempir- 
ical theory of Pollack and Cuzzi [48] was used to derive the 
aerosol phase function as the best fit for the closure between the 
measured ground-based spectral optical data (the aerosol optical 
thickness 7, and the angular sky brightness distribution in the 
almucantar). It should be pointed out that the correct account 
for the aerosol microphysical model (e.g., the shape and size of 
particles) is an important issue of any aerosol retrieval algorithm 
(e.g., see Stamnes ef al. [54]). 

The algorithm has as an input the TOA reflectance at chan- 
nels 1-6 of SeaWiFS (0.412-0.670 um) for the retrieval of the 
aerosol optical thickness over land. Channels 6 and 8 are used 
for the determination of the NDVI and, therefore, the land re- 
flectance. Because of the decreasing surface reflectance with de- 
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Fig. 2. Aerosol optical thickness at the wavelength 0.412 jzm. AERONET 
stations used for the validation are shown by symbols. Black areas correspond 
to clouds. 


creasing wavelength, the channel 1 (0.412 um) gives the best re- 
sults for the aerosol optical thickness retrievals. All eight chan- 
nels are used for retrievals over water. 

The aerosol part of the algorithm was validated by comparing 
with ground-based measurements. The results of satellite and 
ground-based techniques differ not more than by 20% and show 
a sufficient correlation. Details of the validation campaign can 
be found in [56], where also a full description of the algorithm 
is presented. 

The results of the aerosol optical thickness retrieval are pre- 
sented in Fig. 2, which show the spatial distribution of 7, at 
A = 0.412 um at a given time (11:30 UTC) on October 11, 
2001. The clouds are masked as black areas in Fig. 2. The cor- 
responding frequency distributions of the optical thickness is 
shown in Fig. 3. We also show the frequency distribution of the 
retrieved value of a (see Fig. 4). The validation of results using 
AERONET DATA is presented in Fig. 5. 

Let us analyze results obtained. The aerosol optical thickness 
of the scene studied is peaked around the modal value equal to 
0.25 (see Fig. 3), which is close to the value of 0. 27, reported in 
[27] for roughly the same region. The ratio s of the half-width 
of the distribution of Tą, to the modal value of optical thickness 
To is equal approximately to 0.5. The distribution, presented in 
Fig. 3 may be possibly biased due to the presence of subpixel 
clouds (see Fig. 2). This fact and also possible additional illumi- 
nation of a cloudless atmosphere by light, coming from neigh- 
boring clouds (see Fig. 2), may result in the over-estimation of 
the aerosol optical thickness at larger T. 
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Fig.4. Frequency distribution of the Angstrém exponent for the scene studied. 


The regional weather explains the change of Tą across the 
scene in Fig. 2. High values of T4 are present to the south of the 
cloud front system, whereas the small aerosol optical thickness 
is observed to the north behind the large cloud field at the rear 
of the cold front (see Fig. 2). 

Fig. 2 indicates that the aerosol optical thickness is smaller 
behind a cloud system front (over North Sea, Scandinavia, 
Northern Poland and Baltic Sea). This is best explained by 
the cleaning effect of atmosphere by clouds. A significant 
amount of pixels with 7 in the range 0.15...0.25 (green and 
blue color in Fig. 2) is observed over Central Europe, Italy, 
and the Mediterranean. The aerosol optical thickness over Po 
Valley in northern Italy reaches values of 0.3...0.5. This large 
difference (approximately, two times) can be attributed to the 
anthropogenic activity coupled with the local meteorological 
situation. 

We see that our algorithm is capable to quantify the strength 
of sources of atmospheric pollution in Western Europe, using Ta 
as a convenient indicator of the atmospheric turbidity. 

A comparatively large value of the aerosol optical thickness 
is also observed in the left corner at the bottom of Fig. 2. This is, 
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Fig. 5. Retrieved aerosol optical thickness as compared to AERONET data. 
The AERONET locations can be sorted by the value of the optical thickness. 
Then we have the following order: Ispra (the smallest value of Ta), Rome, 
Avignon, Oristano, Venice. 


possibly, due to the simultaneous existence of a water cloud and 
a dust outbreak from Sahara. Note that generally dust aerosols 
and clouds strongly interact, producing a water—mineral aerosol 
mixture. 


In contrast to the optical thickness of cloudy media, the 
aerosol optical thickness can strongly depend on the wave- 
length [20]. Thus, the value of 7, was retrieved not only at 
the wavelength 0.412 um, as shown in Fig. 2, but also for 
several wavelengths in the range 0.412—0.865 um over water. 
For pixels over land, only wavelengths 0.412-0.670 um were 
used. This is due to a strong light reflection by vegetation for 
larger wavelengths, which considerably decreases the retrieval 
accuracy. 

The spectral slope of the aerosol optical thickness was re- 
trieved for each cloudless pixel. As one might expect, it was 
found, that the aerosol optical thickness has different spectral 
slopes a within the scene. The frequency of occurrence of dif- 
ferent Angstr6m exponents & is presented in Fig. 4. Although 
& is biased both over land and sea due to the insufficient in- 
formation on the spectral surface reflectance (especially for the 
channels 0.555 and 0.67 um), it gives an indication for different 
aerosol types. 

A certain fraction of the atmospheric aerosol has nearly neu- 
tral spectral slopes. Similar findings were also reported by Ig- 
natov and Stowe [26] for the case of oceanic aerosols. It means 
that the size of aerosol particles is comparatively large in the 
scene studied here. This is, probably, a consequence of the influ- 
ence of the maritime aerosol or possible cloud formation events 
in the scene under study. The dust aerosol, which originates due 
to the atmosphere—land interactions and due to the transport 
from other regions (e.g., from the Sahara; see bottom of Fig. 2), 
could be another source of small or even negative values of a. 
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Other parts of the scene have values of a, which are typical 
for continental aerosols [20]. The most frequent values of œ are 
found to be in the range [—0.1, 1.1] (see Fig. 4). Note that for 
continental aerosols the value of œ is usually in the range 1.0-1.5 
[2], [28], [51]. Deviations from these numbers obtained for our 
particular case are most probably due to the influence of the 
maritime aerosol and cloud formation processes, as it was men- 
tioned above. We should underline that values of a in the range 
[—0.1; 1.8] were reported for maritime aerosol in [38]. Schmidt 
[51] has obtained values of a in the range [0.1; 0.4] for the clean 
oceanic aerosol. 

Note that operational aerosol retrieval algorithms give values 
of a in the range [0, 3] and even larger [8]. Our results (and those 
in [21], [23], [38], [51], and [52]), however, show that values of 
a larger than 2.0 occur infrequently. This contradiction should 
be resolved in future research. 

Ignatov and Stowe [26] proposed to use the normal distribu- 
tion for the description of the probability distribution functions 
of the Angstr6m exponent over oceans. It follows from Fig. 4 
that the probability distribution function of the Angstrém ex- 
ponent for pixels over land in our case can be also fitted by a 
normal distribution, if we disregard events with large Ta. They 
are most probably due to residual cloud contamination of pixels 
assumed to be clear. 

The frequency distribution of the aerosol optical thickness in 
Fig. 3 is skewed and have long tails for larger Ta. This means 
that they can be approximated by the Gamma or Lognormal 
distributions even for such large scenes as given in Fig. 2. Such 
statistical distributions of aerosol and cloud optical thickness 
have also been reported and used by other research groups [3], 
[12], [25], [26], [46], [47]. 

It is of importance to check the accuracy of our retrieval, 
using ground-based measurements. For this we have used 
data given by the Aerosol Robotic Network (AERONET) 
[11], [19]. The AERONET is the number of identical globally 
distributed sun- and sky-scanning ground-based automated 
radiometers. They provide measurements of aerosol optical 
properties, including aerosol optical thickness on a long term 
basis (up to ten years in some locations). The channels of 
AERONET instruments (0.34, 0.38, 0.44, 0.5, 0.67, 0.87, and 
1.02 um), however, do not coincide with the channels of the 
SeaWiFs (0.402-0.422 um, 0.433-0.453 um, 0.480-0.5 um, 
0.545-0.565 jum, 0.660-0.680 um, 0.745-0.785 um, and 
0.845-0.885 um). So, we made an interpolation between the 
individual spectral points, where the optical thickness is mea- 
sured by five AERONET stations located in Avignon, Ispra, 
Oristano, Rome, and Venice, to have the results for central 
positions of the SeaWiFs bands. Also, we averaged results of 
our retrievals for five pixels surrounding the AERONET sites 
to avoid possible fluctuations in the aerosol spatial distribution. 
Other AERONET stations placed in Europe gave no data due 
to the influence of the cloudiness at their sites (see Fig. 2). 

Results of comparisons for first two channels of the SeaWiFS 
instrument are given in Fig. 5. It follows that the results are 
less accurate for the small optical thickness, which one should 
expect. There is also the deviation for the comparatively large 
optical thickness (close to 0.9; see AERONET results). How- 
ever, we see that the correspondence is generally good for most 


common values of the aerosol optical thickness (around 0.25) 
even accounting for the fact that the AERONET measurements 
are performed at a given point and satellite data refer to much 
broader area. This underlines a comparatively high accuracy of 
our aerosol retrieval technique over land for first two channels 
of the SeaWiFs instrument. We also have compared our satellite 
retrieval results for larger wavelengths (not shown here). It was 
found that the difference of the retrieved optical thicknesses as 
compared to those derived from AERONET generally increases 
with the wavelength. Most probably, this is related to the fact 
that only single spectrum of a green vegetation (measured over 
Germany) is used in our retrieval technique. However, in re- 
ality spectra of vegetation could differ considerably. Note that 
AERONET sites we have selected are located in southern part of 
Europe, where vegetation (and also terrain) differs from that in 
Germany. We would like to address discrepancies found in our 
next publication. The light reflection by vegetation (and also by 
bare soil) is low for first two channels of SeaWiFs. This makes it 
possible to make accurate retrievals of the aerosol optical thick- 
ness over land for these channels. In particular, we note that the 
distribution given in Fig. 3 is representative for the aerosol op- 
tical thickness over Europe. 


HI. CLOUD OPTICAL THICKNESS AND SPHERICAL ALBEDO 


The cloud optical thickness is found using the asymptotic ap- 
proach proposed in [31] and [50]. The retrieval is based on the 
asymptotic solution of the radiative transfer equation at large 
cloud optical thickness. It has the following form [31]: 


where the reflection from underlying surface is neglected. Here, 
H is the cosine of observation angle, jo is the cosine of the 
incidence angle, ¢ is the relative azimuth, Roo is the reflection 
function of a semi-infinite layer, K (p) is the escape function 
[55], and ¢ is the global transmittance of a cloud layer, given by 
the following formula: 


t= (o +0.75r(1 — g)) 1. (2) 


Here, g is the asymmetry parameter [31] and o is the 
parameter, which is approximately equal to 1.07 for cloud 
phase functions [31], [34]. The precise definition of functions 
Roo(Ho, H, $) and K (js) are given in [55]. 

The accuracy of (1) is better than 1% for clouds having the 
optical thickness larger than 10 [31]. Therefore, it provides an 
excellent approximation for the overcast cloudy conditions (see 
the black colored area in upper portion of Fig. 2). King [31] 
validated the algorithm of cloud optical thickness retrieval based 
on (1) and found a good agreement with airborne experiments 
(see also [43]). 

We assume here that clouds are composed of water droplets, 
characterized by the Deirmendjian’s Cloud C1 droplet size dis- 
tribution [9]. As the modified asymptotic theory is accurate only 
at T > 10, large retrieval errors can result for smaller cloud op- 
tical thickness. However, thin clouds often consist of ice crys- 
tals. Also, as optically thin clouds are highly inhomogeneous, 
the application of the standard radiative transfer equation for a 
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homogeneous plane-parallel layer is not justified and may lead 
to large retrieval errors [36], [37]. Generally, radiative proper- 
ties of inhomogeneous clouds are poorly understood. So, this 
study is mostly limited to the case of optically thick clouds. A 
band of a such cloud field is seen in the upper part of Fig. 1. 

The retrieval of the cloud optical thickness is performed at the 
wavelength 0.412 um, where the reflection of light from the veg- 
etation, bare soil, and water can be neglected. However, due to 
the large size of cloud droplets as compared to the wavelength, 
the optical thickness in visible almost does not depend on the 
wavelength. So, results obtained can be safely extrapolated at 
least for the whole visible part of the electromagnetic spectrum. 

The equation for the cloud optical thickness follows from (1) 
and (2) after simple algebraic calculations. In particular, we have 
from (2) 


4tt—o 
o> ce Ss 
where [see (1)] 
t= Rolko, L, P) = Ro, p, $) (4) 


K(uo)K(u) 


Note that the accuracy of the retrieval according to (3) decreases 
as t decreases. The spatial distribution of the cloud optical thick- 
ness derived using (3) is shown in Fig. 6. 

It was assumed that g = 0.85, which is the case for most 
of water clouds. Ro({0, p, 9) was calculated using the code 
described in [40] for the Cloud C1 gamma distribution of water 
droplets [9], having the effective size 6 zm and À = 0.412 um. 
This assumption does not influence the result of the retrieval, 
however. This is because of low sensitivity of the function 
Roo(Ho; H, $) to the size of droplets [34]. Also, we used the 
approximation 


Ku) = 3 


which is valid with the relative error less than 2% at u greater 
than 0.2 [34] for arbitrary phase functions. 

Note in passing that (3) can also be used for the determination 
of the cloud optical thickness from airborne measurements of 
the reflected solar flux [33]: r(u0) = 2 T R(n, po)udp, where 
R(u, po) is the azimuthaly averaged reflection function. Then, 
we have [33]: r(uo) = 1 — tK (mo) for a nonabsorbing wave- 
length and therefore 


(1 + 2p) (5) 


= 1 — r(uo) 
K (wo) l 


So, rather complex calculations of Roo (Ho, H, $) are avoided 
in this case. This is an important point. 

Turning to the analysis of cloud optical thickness distribution 
derived, we note that the cloudiness in the terrestrial atmosphere 
is best described as a persistent and extensive background of op- 
tically thin clouds mixed together with highly variable compo- 
nent of optically thick (and often precipitating) clouds, which 
are very sparse [49]. Namely this case corresponds to the scene 
studied (see Fig. 2), where the band of the optically thick cloud 
field propagates from the North-West to the East in agreement 


(6) 


Latitude 


Longitude 


Cloud Optical Thickness 


Fig. 6. Cloud optical thickness spatial distribution for the scene studied. The 
fractal structure in the middle of the cloud band is clearly seen. 


with other satellite observations (e.g., see the weather map for 
this day at http://www.wetter.zentrale.de/topkarten). 

Note that winds are mostly eastward in midlatitude regions. 
Thus, the maritime aerosol with almost neutral spectral slopes 
are often observed in many areas of the Western Europe. It also 
confirmed by this study. This is also the reason for compara- 
tively large aerosol concentrations, e.g., over Eastern Europe as 
compared to the Western Europe [6]. Larger aerosol concentra- 
tions influence the cloud droplet spectra, precipitation, and con- 
centrations of droplets and crystals in clouds [6], [44]. 

The study of the cloud field shows that the modal cloud op- 
tical thickness in the brightest part of the cloud band in Fig. 6 is 
around 30 (see Fig. 7). 

Note that (4) can be used for the calculation not only the 
global cloud transmittance but also the cloud spherical albedo 
r = 1 — t, which is an important parameter for climate studies. 
It gives a hemispherical reflectance of a cloud layer. Then the 
additional information on the value of g is not required to run 
the retrieval procedure. This is an important point. 

We have derived the spatial spherical albedo distribution over 
the cloudy portions of the image studied. This distribution (not 
shown here) is similar to that of the cloud optical thickness 
given in Fig. 6. The frequency distribution of the cloud spher- 
ical albedo in the scene (see Fig. 8) differs considerably from 
the correspondent frequency distribution for the cloud optical 
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Fig. 7. Cloud optical thickness frequency distribution. 
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Fig. 8. Cloud spherical albedo frequency distribution for the scene studied. 
thickness shown in Fig. 7. The maximum of the spherical albedo 
distribution is located around 0.8 for the case under study. 

The distribution itself can be described by the well-known 
beta distribution, which corresponds to similar findings in [35]. 
However, their results are based on pixel brightness studies for 
a given solar angle and observation geometry. Our results, pre- 
sented in Fig. 8 correspond to the angle-independent spherical 
albedo. Note that the betta distribution usually appears if there 
is a boundary for the physical value studied. In our case, for ex- 
ample, cases with r > 1 are not possible due to physical reasons. 


IV. CONCLUSION 


The atmospheric optical thickness over land and water sur- 
faces has been retrieved, using newly developed aerosol-cloud 
retrieval algorithm. 

Note that our algorithm is not capable to deal with partly 
cloudy pixels. The cloudy part of the joint algorithm works only 
for clouds having optical thickness larger than approximately 5. 
The aerosol part of the technique can be applied for studies of 


aerosol optical thickness distributions both over land and ocean. 
However, land surfaces with large reflectance (e.g., snow fields 
and deserts) should be avoided. Then the aerosol signal is too 
small as compared to the surface contribution. This does not 
allow to get meaningful results using SeaWiFS imagery. 

The most probable cloud optical thickness observed for the 
extended cloud field in Fig. 2 is around 30 (see Fig. 7). Note also 
that one can observe a fractal structure in the middle of the cloud 
band in Fig. 6. The most probable value of the aerosol optical 
thickness is equal to 0.25 (see Fig. 3). Note also the large vari- 
ability of 7, and regional differences in the scene (see Fig. 2). 

The results reported in this paper have a high accuracy for 
large (T > 10) optical thickness. Generally speaking, lookup 
tables should be used in the case of thin clouds. However, the 
question arises if the exact radiative transfer theory for plane- 
parallel layers can be applied in the case of thin clouds at all. 
This obstacle is due to the pronounced stochastic nature of thin 
clouds. The problem should be checked and clarified in future 
research (see [3]-[5], [36], [37], and [39]). We have also re- 
trieved the cloud optical thickness outside the extended cloud 
field in the upper portion of Fig. 6. However, results obtained 
are most probably biased due to errors in (1) for the smaller 
values of the optical thickness. 

The data presented here (e.g., see Figs. 2 and 6) can be 
used for the comparison of different aerosol/cloud retrieval 
algorithms, existing worldwide. Such studies are needed to 
overcome the current complexity, associated with the theory 
and practice of the cloud and aerosol remote sensing in the 
terrestrial atmosphere. 

Note that application of our algorithm reveals some important 
issues to address in future. First of all we note that the aerosol 
optical thickness tends to be uniform outside industrial areas and 
special events like fires, outbreaks from deserts and volcanic ac- 
tivity. Therefore, highly variable aerosol optical thickness dis- 
tribution shown in Fig. 2 is, probably, due to the cloud formation 
influence. It means that further study of the cloud threshold in 
the aerosol retrieval scheme is needed. Further, the spatial dis- 
tribution of the Angstrém exponent (not shown here) differs sig- 
nificantly from the average data for continental areas (a = 1.3 
[28]). This is, probably, due to the influence of maritime aerosol 
and cloud formation processes taking place in a given scene as 
well. On the other hand, it could be also partially due to the sur- 
face variability, which is taking into account only approximately 
in our algorithm. Our surface reflection model should be further 
improved, therefore. 

Due to the limited spectral range of the SeaWiFS instrument it 
is not possible to retrieve other cloud parameters such as cloud 
top height, cloud liquid water path, cloud droplet effective ra- 
dius, cloud geometrical thickness, cloud thermodynamic phase 
(and similar parameters for the atmospheric aerosol). We also 
have worked out the extension of our algorithm to treat these 
more complex retrievals, using spectral distribution of backscat- 
tered light for cloudy scenes in visible and near-infrared spectral 
ranges (e.g., see [34]). 

Our next task is to apply our retrieval scheme to aerosol— 
cloud media in different geographical locations for different 
seasons. This will enable us to gain an important statistical 
information on cloud—aerosol interactions and properties. 
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In conclusion, we note, that extensive satellite retrievals of at- 
mospheric turbidity given by 7 (and also other aerosol and cloud 
micro- and macrophysical parameters, e.g., particle size dis- 
tributions and shapes, cloud-top-heights and cloud geometrical 
thicknesses, cloud thermodynamical state, etc.) provide the sta- 
tistical data on the required temporal and regional variability of 
suspended matter in the atmosphere of our planet. This informa- 
tion along with the data on the concentrations and variability of 
trace gases [7] is of a paramount importance for modern atmo- 
spheric chemistry and physics as well as climate change studies 
[24]. 
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